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Abstract 

We present the solution of the problem of the !/□,□ — 0, asymptotic 
terms discovered in the one-loop form factors of the gravitational effective 
action. Owing to certain constraints among their coefficients, which we es- 
tablish, these terms cancel in the vacuum stress tensor and do not violate the 
asymptotic flatness of the expectation value of the metric. They reappear, 
however, in the Riemann tensor of this metric and stand for a new effect: a 
radiation of gravitational waves induced by the vacuum stress. This coher- 
ent radiation caused by the backreaction adds to the noncoherent radiation 
caused by the pair creation in the case where the initial state provides the 
vacuum stress tensor with a quadrupole moment. 
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1 Introduction 



We consider the expectation-value equations for the gravitational field in an 
in-state [} The model-independent, or phenomeno logical approach [4-7] makes it 
possible to write down the general form of these equations in terms of the form factors 
in the vacuum action. The form factors are to be calculated from a given dynamical 
model. However, for obtaining predictions of various models, the expect at ion- value 
equations should be first analysed with arbitrary form factors in order to relate the 
properties of the form factors to the important properties of the solution [4]. 

This analysis has thus far been limited to the behaviour of the vacuum stress 
tensor T^"^ at infinity [6,7]. It was shown [6] that the requirement of asymptotic 
flatness of the solution imposes restrictions on the behaviours of the form factors 
at small values of their arguments. Namely, these behaviours in one (each) of the 
arguments with the others fixed should be t/7(0) log(— □) + 0(1), □ —>■ —0. The 
coefficients w{0) of the log(— □), □ —>■ —0 behaviours (the spectral weights at zero 
spectral mass) determine the rate of the vacuum radiation through the future null 
infinity (X+). 

However, at each given order in the curvature, only certain combinations of the 
form factors must behave in this way. The condition of finiteness of T^^^ at X"*" which 
brought to the result above leaves some arbitrariness in the asymptotic behaviours 
of the individual form factors in the basis decomposition of the action. It is this 
arbitrariness that allows for the existence of the effect discussed in the present paper. 

The present paper deals with the problem which appeared when the field-theoretic 
form factors were calculated. It is worth noting that, as emphasized in [6], the loop 
expansion of field theory has a domain of validity. It is near null infinity where the 

^ The gravitational collapse problem was first considered in this setting in refs. [1-3]. 
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results of this expansion are valid and can be used to calculate the energy of the 
vacuum radiation. The ultraviolet divergent terms which appear in the expectation- 
value equations when T^*^ is expanded in loops are local and vanish at infinity; 
only nonlocal terms survive, and these are unambiguous. Technically, the loop ex- 
panded form factors are reliable in the limit of small □ arguments up to terms 

O(n0),n ^ -0. 

The field-theoretic form factors have been calculated in the one-loop approxima- 
tion for a generic quantum field model in refs. [8-15]. Their asymptotic behaviours 
at □ — > — up to terms 0(n°), which are of interest for the above-mentioned reason, 
are presented in ref. [15], and these behaviours offer a problem. While the second- 
order form factors behave as expected: they are w{0) log(— □) + 0(1) for □ — 0, 
the third-order form factors contain also the asymptotic terms l/D, □ — >• —0 which 
apparently violate the asymptotic flatness of the solution [6,15]. 

Since the third-order form factors are functions of three □ arguments, their 
behaviours in one of the arguments with the two others fixed cannot be predicted 
on dimensional grounds. The l/D asymptotic terms in the form factors appear as a 
result of an explicit calculation of loops [12,15]. In the present paper we propose an 
explanation of this result as well as of the following remarkable fact [12] which one 
can establish by a direct inspection of the expressions in [15]. The inspection shows 
that the alarming I/O terms appear only in the curvature invariants containing 
the gravitational field strength and act selectively only on the Ricci curvature. The 
matter field strengths contained in the commutator and potential curvatures (see 
[15]) remain unaffected by these terms. Thus the presence of the l/D terms breaks 
the democracy of massless vacuum particles; gravitons appear to be distinguished. 

Since the same vacuum action describes also the transition amplitudes between 
in- and out- states [8] , the problem of the I/O terms appears also in scattering theory 
where these terms either signal an infrared divergence of the on-shell amplitudes 
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with gravitons or, in the favourable case, stand for some inelastic process. In the 
language of expectation values, the former case corresponds to a breakdown of the 
asymptotic flatness. That the situation is not hopeless is seen from the fact that the 
form factors are not quite the vertices; they are coefficients of the curvatures rather 
than the field disturbances. As compared to the field disturbances, the curvatures 
contain extra derivatives which, in the favourable case, may cancel the !/□ factors 
in the on-shell amplitudes. 

Below we present the solution of the problem as it appears in expectation- value 
theory. As mentioned above, only certain combinations of the form factors should 
behave hke log(— □), □ ^ — to ensure finiteness of T^^^ at J"*". We show that the 
!/□ terms precisely cancel in these combinations leaving indeed the log(— □) be- 
haviour as the leading one. The cancellation occurs owing to certain constraints 
among the coefficients of the !/□ terms which we establish by an analysis of the 
asymptotic flatness and next check with the explicit expressions in [15]. The fulfll- 
ment of these constraints is by itself a powerful check on the results in [15] apart 
from the checks that have already been carried out in [12]. Thus we prove that the 
!/□, □ — > —0 terms discovered in [12,15] do not violate the asymptotic flatness of 
the solution of the expectation-value equations. The proof is given for a generic 
quantum fleld model for which the results in [15] are obtained and which is charac- 
terized by a set of fleld strengths consisting of the Ricci, commutator, and potential 
curvatures. Such a general proof is possible owing to the above-mentioned fact that 
the 1/n operators in the asymptotic expressions for the form factors act only on the 
Ricci curvature. Therefore, for the consideration of the leading asymptotic terms of 
the equations at X+, one does not need to know the variational derivatives of the 
commutator and potential curvatures with respect to the metric. 

Next, we reveal the signiflcance of the l/D, □ —0 terms in the vacuum form 
factors. These terms vanish in the energy-momentum tensor but , as we show, they 
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reappear in the Riemann tensor of the solution and determine its leading 0(l/r) 
behaviour at X"*". The coefficient of this behaviour is known to give the energy of 
the outgoing gravitational waves. Thus the !/□, □ ^ —0 terms of the vacuum form 
factors discovered in [12,15] stand for a new effect: a generation of the gravitational 
waves from the vacuum. This is an effect of the backreaction of the vacuum stress 
on the metric. All massless particles including gravitons [16] contribute to Tj^^ 
and are radiated through the future null infinity by the quantum mechanism of 
pair creation. The energy of this component of radiation is determined by the 
log(— □), □ — > — terms of the vacuum form factors. On top of this, T^^ as a 
whole acts as a source of the gravitational field and causes a secondary radiation of 
gravitons. This component of radiation has the shape of a classical wave but with 
a quantum amplitude, and its energy is determined by the !/□, □ — > —0 terms of 
the vacuum form factors. 

An important difference between the two cases is that the gravitational- wave 
component will be nonvanishing only if the initial state has a sufficient asymmetry 
to provide the vacuum stress with a quadrupole moment whereas the contribution 
of gravitons in T^ac\i+ present even in a spherically symmetric in-state because 
the out-states in which these gravitons appear at I'^ are squeezed vacuum states 
rather than coherent states (see [17,18] and references therein). 

The paper is organized as follows. In sec. 2 we briefly review the structure of 
the expectation- value equations and present a new expression for the solution of the 
Bianchi identities to second order in the curvature. This expression simplifles ob- 
taining the news functions of the gravitational waves. Sec. 3 contains an analysis of 
contributions to at X+. This analysis is not complete but sufficient for obtaining 
the asymptotic flatness constraints and for calculating T^'^|x+ to second order in the 
curvature. We point out an important distinction of this calculation from the one 
in two dimensions [3] , and also make some flrst step in giving the expectation- value 
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equations a closed form. Sec. 4 contains the derivation of the asymptotic flatness 
constraints and the proof that they are satisfied by the one-loop form factors. In sec. 
5 we propose a new method for calculating the energy of the gravitational waves 
and obtain the contribution of the third-order vacuum form factors to the news 
functions. Sec. 6 completes the calculation of the vacuum news functions in the 
lowest nonvanishing approximation. Appendix A contains reference equations per- 
taining to the behaviour of the asymptotically flat metric at null infinity. Appendix 
B summarizes the properties of the retarded Green function used in the text. 



2 The expectation-value equations in an in-state 



In the framework of quantum field theory one starts with the assumption that 
there exists a quantum state such that the expectation value of the metric in this 
state is an asymptotically fiat gravitational field. Under this assumption one goes 
to the past null infinity (X^) of the spacetime with the expectation value of the 
metric and, for all massless fields, builds the Fock space of states (the in-states) 
based on the standard in- vacuum [19]. The assumed state belongs to this space 0. 
The choice of the state determines the initial data at for the field's expectation 
values and, generally, affects also the form of their dynamical equations since these 
equations are nonlocal. There exists an action which produces the expect at ion- value 
equations although the procedure by which it does so is not the least-action principle 
(see below). 

^After the choice of the state has been made, one checks the original assumption. Since the 
causaUty relationship is set by the expectation value of the metric, this is a self-consistent problem 
even in its original setting. That the consistency check is nontrivial is seen, for example, from 
the fact that with the massive quantum fields one generally arrives at a contradiction with the 
asymptotic flatness. 
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It makes sense to choose the initial state in which the matter quanta form some 
heavy classically behaved source of the gravitational field, and gravitons are in a 
coherent state so that, for the mean metric, they form, generally, a classical incoming 
gravitational wave. Here we consider the case where such a wave is absent. The 
action for the expectation value of the metric in such a state can be taken as the 
sum 

^ ^vac ~l~ ^source (^■-^) 

where Ssource is the action of a source which moves along a classical trajectory in 
the mean metric, and S^ac is the action for the gravitational field in the in- vacuum 
state. 

The action S^ac is to be calculated from a given quantum field model. Within 
certain approximations (which are not completely unsatisfactory, see above) this 
calculation is feasible and, in the one-loop approximation, it can be done for a 
generic field model [8-15]. However, since the gravitational interaction is universal, 
all particles existing in nature contribute to S^ac, a-nd, at higher loop orders, all 
details of their interactions matter. Therefore, if S^ac is to be ultimately calculated 
from a model (of fields or strings or whatever) , then this should be the Model and 
the Calculation. 

In the phenomenological approach of refs. [4-7], the action (2.1) is viewed as an 
effective action (in the loose sense) for the observable field which should be a part 
of predictions of any fundamental dynamical theory. Irrespectively of the nature 
of this theory, one assumes the existence of a functional, the action S^,ac, which 
describes the elastic properties of real vacuum i.e. its response to the introduction 
of a gravitationally charged source. For Syac one writes down the most general 
expansion in terms of nonlocal invariants of A'"th order in the curvature. One has to 
go explicitly to = 3 because third order in the curvature in the action corresponds 
to second order in the equations, and it has been shown that, at first order in the 
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curvature, the flux of vacuum energy tlirougli X+ is pure quantum noise [7]. Tlie 
full bases of nonlocal invariants of second order and third order are built in [5] for a 
set of field strengths consisting of the Ricci, commutator and potential curvatures: 

3? = {i^^,,7^^,,p} (2.2) 

( for the definitions see [15]). The explanation for the absence of the Riemann tensor 
from the basis invariants can be found in [9,4,5] but it makes sense to repeat it here 
since the respective equations will be of use below. 

By differentiating and contracting the Bianchi identities, one obtains the equa- 
tion 

- R^^^^R^"'""' (2.3) 

which can be solved iteratively with respect to the Riemann tensor. In this equation 
the Ricci tensor plays the role of a source, and the solution is fixed by the initial 
data for the gravitational field at I~ . The solution with zero data (no incoming 
gravitational wave) corresponds to the in-vacuum state and is expressed in terms of 
the retarded Green function (see appendix B): 

= l(^4V[^V<"i?''l'^> + 0[R^]). (2.4a) 

Here and below, !/□ stands for the retarded Green function, and both types of 
brackets [] and <> denote antisymmetrization of the respective indices. Since the 
Riemann tensor is expressed in this way through the Ricci tensor, the nonlocal 
invariants with the Riemann tensor in the vacuum action are redundant. 

Below we shall confine ourselves to the case where the flux components of the 
Ricci tensor at vanish |^. In this case the derivatives in (2.4a) can be made 
^ For the expectation value of the metric this condition restricts the choice of the classical source 
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external by commuting them with the Green function !/□. To second order in the 
Ricci curvature the solution is then of the form 



(2.4b) 

+ 



+ 4(-V['^V^i?''l<") {-R'f) + 4V['^([-V<"i?^>^ 



□ 



(cf. the result in [5,12]). The advantage of making the derivatives external is in the 
appearance of the terms in the curly brackets in (2.4b) which have no overall l/D 
factor. These terms do not contribute to the leading asymptotic behaviour of the 
Riemann tensor at X+. On the other hand, the terms which have the overall l/D 
factor have also the overall derivatives. This facilitates solving the equation for the 
news functions of the gravitational waves ( see sec. 5). 

To third order in the curvature the vacuum action is of the form ^ 



S,ac = Sil) + Si2) + S{3) + 0[^' 



S{1) 



167r 



3(2) 



1=1 



(2.5) 
(2.6) 

(2.7) 



which can be taken, e.g., as having a compact spatial support. Then, since the energy-momentum 
tensor of the in- vacuum has no incoming fluxes at , the same wifl be true also of i?^'' (see eq. 
(2.18) below). In the general case, the flux components of R'^'^ at cancel in the combination 
appearing in (2.4a) but the derivatives cannot be commuted with !/□ for otherwise 

the action of the retarded Green function will become ill-defined. 

''We follow the notation of ref. [15] but change the overall sign of the action as appropriate 
for the lorentzian signature of the metric. We use the signature ( — h ++) and the conventions 
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1 r 

"^(3) = ^^77^ / c?2;^'/'tr^r,(-ni,-n2,-n3)3?i5R23?3(i), (2.8) 
2(47r) 

□ = ^''"V^V. (2.9) 

where 3?i3?2(0 with i = 1 to 5, and ^i^2^z{i) with i = 1 to 29 are the quadratic 
and cubic basis invariants hsted in [15]. This hst is reduced in comparison to the full 
list in [5] because, in the trace of the heat kernel and hence in the one-loop vacuum 
action, the invariants linear in the commutator curvature all but one prove to be 
absent [12,13]. The only one that is present is number 13 in the list of rcf. [15] which 
we use here. In low-dimensional manifolds there exist hidden constraints between 
nonlocal invariants, reducing the basis. In four dimensions, the second-order basis is 
unconstrained, and the only constraint which exists among the third-order invariants 
boils down to the condition that the completely symmetric part of the form factor 

vanishes identically [5,12]. In the field-theoretic form factors of refs. [12,15] this 
condition is explicitly implemented. 

The commutator and potential curvatures are functions of the metric and matter 
fields, different in different models, but in any case their contribution to the purely 
gravitational sector of the action boils down to a modification of the form factors of 
the basis invariants with the Ricci tensor only. There are only two such in S{2) and 
ten in 5'(3). Below, when referring to the purely gravitational form factors, we shall 
assume that this reduction has already been made. The full set of invariants for 
gravity and matter is considered here because it is important that the maintenance 
of asymptotic flatness be proved for the one-loop action in full generality. Apart 
from this, our main concern in the discussion below is the vacuum action for the 
metric. 

The functions % and Fj in (2.7) and (2.8) are the second-order and third-order 
form factors. The principal assumption about these and higher-order form factors 
made in the axiomatic approach is their analyticity which allows one to put them in 
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the spectral forms [4] . For example, the spectral form used in [6] for the lowest-order 
form factors is 



■)<-□) = cvr / +'j:t}iin+,r(i^] i(p') (2.10) 



where w{im?) is the spectral weight 



wirn?) 



1 



7(— m — iO) — 7(— m + iO) 



[2.11] 



27ri 

/i^ > is an arbitrary parameter on which 7(— □) does not depend, and n is the 
degree of growth of 7(— □) at large □ which will presumably be fixed or bounded by 
the requirement of regularity of the solution (see [1-3]). The requirement of asymp- 
totic flatness of the solution imposes restrictions only on the small- □ behaviours of 
the form factors [6,7]. In the small-D limit, eq. (2.10) reduces to the simple spectral 
form since the terms modifying this form for n > vanish in this limit [6]. On the 
other hand, this form should be generalized to allow for the behaviour l/D, |n| — > 
of 7(— □). It remained unnoticed in paper [6] that the derivation in this paper for 
the two purely gravitational second-order form factors brings in fact to the following 
general result: 

7i(-n) = -^-^/;i(0)log(-n)+O(l), -D^O (2.12) 
72(-n) = ^-W2(0)log(-n) + 0(l), -D^O (2.13) 

in which there appears an arbitrary constant a 0. Only the combination 

7i(-n) + 272(-n) = -(wM + 2^i;2(0)) log(-n) + 0(1), -□ ^ (2.14) 

a, — □ ^0 which had 



^Eq. (39) of ref. [6] admits one more solution: □ 7i(— □) + 372(— □) 
been overlooked. 
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should behave hke log(— □) by the analysis in [6]. As will be seen below, a similar 
situation takes place for the higher-order form factors. 

In the form factors 7^ calculated from field theory, a = [9,15] since, by dimen- 
sion, the terms l/D cannot appear in the loop expansion of the second-order form 
factors. They appear, however, already in the third-order form factors Fj [15]. For 
the generalized spectral forms of these form factors see [10,12]. 

With the form factors in the spectral forms, the only nonlocal structure that 
remains in S^ac is the inverse operator l/(m^ — □). This simplifies the procedure 
of obtaining the expectation-value equations. When the action S^ac is varied, the 
inverse operators are regarded as obeying the variational rule 

^ 11,1 , , 

-.^^-^ ' (2-15) 



— □ — □ — □ 

and, after the variation has been completed, all inverse operators are replaced by 
the retarded Green functions If 



(2.16) 



is to denote the result of this procedure, and a further notation is introduced to 
separate the classical term of Syac '■ 

2 SS^if 



rp^U _ 

vac - §g 



' ''^'^ (2.17) 



then the expectation- value equations corresponding to the action (2.1) are of the 
form 

j^,u _ 1^,.^ ^ 8^^^,. ^ ^2.18) 

where 

rpiiu _ 2 6Ssource 1q^ 
source ^ /I X ' y^'i-J) 



®For the derivation of this procedure in QFT see [8] and references therein. In the phcnomeno- 
logical approach this set of rules is taken for granted [4] . 
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and T^^^ in (2.17) can be interpreted as the energy-momentum tensor of the in- 
vacuum. These equations are to be solved with zero initial data for the gravitational 
field at X~. One arrives at a Cauchy problem [2,3] for nonlocal equations with the 
retarded kernels which are to be integrated from X~ to the future until the solution 
hits a singularity if there remains one. One hopes that it doesn't. 

We do not consider the more general initial data, with a gravitational wave at 
because in this case the action should also be calculated more generally. Specifically, 
the solution of eq. (2.3) can no more be taken in the form (2.4). 

3 The structure of T^f^^ at X+ 



It is natural to begin the study of the expectation-value equations with the 
behaviour of T^'^ at null infinity since the nonlocal terms of the equations should 
be responsible for the effect of the vacuum radiation. 

At the future null infinity one has the Bondi-Sachs equation [20,21] (see also 
appendix A) which is an exact consequence of the expectation- value equations : 



dM{u) 
du 



1 

47r 



'd_ 
du 



du 



Co 



1 



1+ 



(3.1) 



(Vm)2 = 0, (Vm, Vr) 
(Vv)2 = 0,{Vu,Vv] 



1+ 



-1, 

: -2 



(3.2) 
(3.3) 



Here u is the retarded time along X+ with the natural normalization in (3.2), the 
integrals are over the 2-sphere S (normalized to have the area An) at which the 
null congruence u =const. crosses X+,r is the luminosity distance along the rays of 
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this congruence, M{u) is the Bondi mass, and dCi/du, dC2/du are the Bondi-Sachs 
news functions of the gravitational waves. 

Eq. (3.1) is the conservation law missing in the theory of quantum fields on a 
fixed gravitational background. In the collapse problem, this is the backreaction 
equation relating "the changing mass of the black hole" with the energy of the 
quanta radiated by this black hole. In full quantum theory, both the "black hole" 
and the quantum fields "on its background" evolve from one and the same initial 
state, and one is able to answer the question where does the black-hole radiation 
take its energy from. It takes it ultimately from the energy of the collapsing source 
'^source which equals the ADM mass of the expectation value of the metric and serves 
as an initial datum M(— oo) for eq. (3.1). 

The last term of eq. (3.1) is the flux of the vacuum energy through . For it 
to be finite, the fiux component of T^^^ should decrease at X"*" like 1/r^. This is a 
necessary condition of asymptotic flatness. Below, terms 0(l/r^) in T^^^ will be 
referred to as vanishing at X"*". 

When computing T^^^{x) at X+ from the action (2.5), all terms in which the 
Ricci curvature appears at the observation point x can be discarded because Rf^i, 
decreases at least hke 1/r^ and will always be multiplied by a decreasing function. 
Thus, the covariant derivatives V which appear in the basis invariants 3?3?3?(i), etc. 
need not be varied because the contributions of their variations to T^*^ contain the 
curvature at the observation point. For the same reason, in the expression 

5KI = ^9''^{w^W^5g,a + V.V"%a - V^V.^^^^^^^^ - + 

+ 2R;^^^6g^f,) + ^R^'-Sg.a - Ig'^'R^jQap , (3.4) 

the terms with the Ricci tensor can be discarded but the term with the Riemann 

tensor cannot since the Riemann tensor has components decreasing like 1/r. In (3.4), 
the expression for 6R2 has been brought by commutations to the form used below. 
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Only the variations of the Ricci tensors and the variations of the form factors 
in Svac can give nonvanishing contributions to T^'^ at . It is easy to see that the 
variation of a form factor 

r(-ni,-n2,...,-niv) (3.5) 

in the argument Dp can contribute at only if, in this argument, the form factor 
behaves like l/Dp, Dp —0. Thus, assuming a = in eqs. (2.12), (2.13), as is the 
case in the field theoretic form factors, one can calculate the variations of 7j(— □) 
in the action (2.7) by using the spectral form (2.10) in which, moreover, the terms 
appearing at n > can be disregarded when the observation point tends to X+. One 
finds 



oo 

J J J TTXi I— I TTh 



□ 



Since, by the result in [6] , 



{m? - D^et) 



OC r ^{x) exp^— |const.|my^r(a;)j (1 + C), 
e)^0,r(x) ^oo,x^X+, 



(3.6) 



(3.7) 



and 



(3.8) 



w(0) = finite, 

one concludes that the contribution of 5^{—U) to T^^^^ at X+ is 0(l/r^). This result 
applies to all form factor of the form (3.5) provided that their behaviours in indi- 
vidual arguments are 0(log(— □)), □ — > —0. On the other hand, with the behaviour 



r(-ni, -02, -Da, ...) = —Fin^, Dg, ...)(i + O), O ^ o, Di 



the variation of Di in the term 

' dxg^'^T{-U^, -02, -Ds, ...^lU^Uz ... 

16 



I' 



(3.9) 



(3.10) 



of the action is asymptotically of the form 

- J dxg'/\^^i)6a^[F{a,, Dg, ...)^,^, . . .) (3.11) 

and gives a nonvanishing contribution to T^^^ at X"*" proportional to 1/r^. 

The fact that, in four dimensions, the variations of the second-order form factors 
do not contribute to the energy flux at infinity makes an important distinction of 
this case from the case in two dimensions where 7(— □) oc 1/n and the relevant 
contribution comes from 6'~f{—0) [3]. The variations of the third-order form factors 
Fj in (2.8) can already contribute to T^"^ at X+ owing to the 1/n, □ —0 terms 
discovered in [15] but, to second order in the curvature (in the equations), only the 
Ricci tensors in the action S'(3) need to be varied. 

For all terms of the vacuum action, in which the form factors are of the form 
(3.5)0, it is useful to introduce a quantity, the generalized current, defined by varying 
the action with respect to the Ricci tensor only. The sum of such terms in S^ac can 
be represented in the form 

S^ac=^^,T.j:SriN), (3.12) 

S,{N) = J dxg'/^P,{Vi). . .Pjv(V,v)r(-ni, -n^)RXx,). . .R:{x^) (3.13) 

where (3.13) is a contribution of the A^th order basis invariant number " i " , and 

the sum in (3.12) extends over both i and A^. In (3.13),-R; are the Ricci tensors 

with mixed indices, and the polynomials in covariant derivatives P„(V„) which are 

generally present in the tensor basis [15,5] act on the respective R.{xn) after the 

action of the operator arguments □„ of F. In this representation all operators □„ in 

F are uniformly defined as applied to a mixed second-rank tensor, and the advantage 

We do not consider here the more general form factors r(— Di, — Dat, —□1+2, — ni_|_3, ...) 
in which the operator arguments □n+m act on products of two curvatures since such form factors 
appear in S^jac only beginning with = 4 [4,5]. 



17 



of taking it mixed is in the absence of g^'^ factors contracting indices in (3.13) 
which otherwise would need to be varied. These factors are generally present in the 
polynomials P(V) but, in any case, their variations ( as well as the variation of g^^'^ 
in the measure) do not contribute to Tj^^ at T"*". The order of operations in (3.13) is 
different from the one in (2.8) (cf. the exphcit expressions in [15]) but, at every order 
in the curvature, the action can be brought to the form (3.13) by commutations. 

Since the argument of F is the first operator acting on R.{xn) in (3.13), it 
will be the last in the variational derivative of (3.13) with respect to R,{xn)- The 
generalized current I^{^,x) is then defined by the relation 



SnS^ac = / dxg'/'lj^i- n,x)SR;{!c)= J dxg"Hi:{- n,x)SR;{x) 



{Any 



(3.14) 

where the notation points out that only the Ricci tensors in (3.13) are varied, and 
the argument □ of J^(— □, x) is the argument of the form factor F in (3.13) that acts 
on the varied Ricci tensor. The I^{$,,x) is a tensor function of the spacetime point 
X and a function of a parameter ^ which in eq. (3.14) gets replaced by the operator 
— □. This operator next acts in either of the two ways pointed out in (3.14). 

Given the action of the form (3.12), it is easy to calculate /^(^, x) to each given 
order in the curvature. For the action (2.5) we have 

/^(e, x) = /^^(e, x) + x) + 0[^] (3.15) 

where /2^(C, x) and /3^(^, x) are the contributions of S{2) and 5'(3) respectively, and 

I^%i,x) = 27i(0^^(^) + 272(0^^^(^) , (3.16) 



+ 



Si^r.oiC, -a^)R';^{x^)R^X2) + 
Vi"V2^F22(^, -a^)R{xi)R{x2) 



Xl=X2=X 



(3.17) 



V^T22(-ni, -□2,0^(^(2^1)^(2^2) 



+ . . . 



X1=X2=X 
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where only the purely gravitational terms are written down. We do not present the 
latter expression in full but exemplify it with the contributions of two third-order 
invariants, number 10 and number 22 (see [15]). Eq. (3.17) illustrates the general 
structure of the current I^{^,x). 

A remarkable property of /^(^, x) is that the variations of both the Ricci tensors 
and the form factors in (3.12) are expressed entirely through this quantity. Indeed, 
in addition to (3.14), we have 

^rS.ac = ^ / dxg'/'Sn'^±^[i;{-n',x") - I;{-D",x")]r^^{x') (3.18) 

where the notation 5r points out that only the form factors F in (3.13) are varied, 
and we have used the general formula for a variation of an operator function [12,14]: 

/ dxg'/'A {6f{n)) B= f dxg'/'Sn^l^^^^^^^^AB . (3.19) 

It is understood that (or □ b) is the operator □ acting on A (or B) , and similarly 
in (3.18) acts on x', and □" on x" with subsequently setting x' — x" — x. The 
operators and n^, and similarly 50' and in (3.18), do not commute and act 
in the indicated order. The identity 

J dxg^/^SD" - 5n')F{x',x") = J dxg^/\-5 log g^/^){n" - n')F{x',x") (3.20) 

(with an arbitrary two-point tensor F{x',x") contracting into a scalar at x' — x") 
serves to check that varying the left-hand side of the equality 



J dxg'/^f{ns)AB = J dxg'/'f{aA)AB (3.21) 



with the aid of eq. (3.19) gives the same result as varying its right-hand side. 

Although the generahzed current I^{^,x), as calculated from the action (2.5), 

is given in the form of an expansion, it enters the expectation-value equations as a 
single whole and determines the vacuum stress at null infinity. Indeed, since (3.14) 
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and (3.18) are the only contributions surviving in (2.16) when the observation point 
tends to X+,T^^ at X+ is obtained in a closed form. As seen from (3.16), (3.17), 
the behaviour of /^(^,a;) in ^ includes all the behaviours of the form factors in 
individual arguments. The I^{^,x) may, therefore, be of a significance in axiomatic 
theory. For this current as a function of its parameter argument one must postulate 
the existence of a spectral form similar to (2.10) but allowing for the spectral weight 
to have a 5(m^) singularity at — 0. Assuming for simphcity n = in (2.10) (the 
modification concerning the large □ is irrelevant to the present discussion) we set 



dm? 



rn? + ^ 

— e 

and then eqs. (3.14) and (3.18) take the form 

1 



^/;^(m^x) , 



5rS, 



T'-'vac 



2(47r)2 

1 °° 
--jdxg^'^jdm 



m? — □ 



(3.22) 



<(m^a;) , (3.23) 



2(47r) 



5U- 



—Ki^) 



□ 



w^lm^, x) 



(3.24) 



Along with expression (3.4) and an easily derivable expression for 5n in (3.24) they 



determine T^^^ at X"*". 



4 The asymptotic flatness constraints 



Even with the 1/^, ^ — > behaviour of x), the contribution to T^^ coming 
from (3.24) is 0(l/r^) at X+. Also the contributions coming from the term with the 
Riemann tensor and the term with the □ operator in (3.4) are 0(l/r^). We have, 
therefore, from (3.14) and (3.4) 

I 1 I J 
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I 1 

- g^'V^VJ^^^i- n.et, x))+ 0{r-\x)), X ^ J+ 

/"^(e,x) = (7("^(x)/f)(e,x). (4.2) 

The terms of order r^^(x) in (4.1) are the ones to be retained but it is not our 
purpose here to calculate T^"^ 0. The point is that the remaining terms in (4.1) 
should also behave like 0(r^^(x)) for the asymptotic flatness to be maintained. 
This behaviour should , moreover, hold at every order in the curvature because even 
small disturbances of the metric can violate the asymptotic flatness. This does not 
mean, however, that the function 

I 1 

0,,,,x) (4.3) 

should behave like 0(r~^(a;)); it suffices that 

I 1 

V„r^(- □,et,x ) = 0{r-\x)) , x^I+. (4.4) 

Owing to this fact, the 1/.^, ^ behaviour of x) is not completely ruled out 

but condition (4.4) imposes a constraint on the coefficient of this behaviour. 

For implementing this constraint the derivative Vq, in (4.4) should be commuted 
with the operator Dret : 

V^r^i- a,et,x) = I dniH ^ ^ VaW^im^x) + (4.5) 

Here all operators act to the right on x, and 

w^^m", X) = g^^\x) w]:){m\ x) . (4.6) 
® For this calculation to lowest order in the curvature see [6,7]. The cancellation of the l/D 
terms being established in the present paper, this calculation can now be done to second order in 
the curvature on the basis of the results in [15]. 
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Both terms in (4.5) are generally 0{r ^{x)) but the commutator term contains an 
extra power of the curvature. If we denote 



, x) = V„w"^(m2, x) + [Va, □] , ^ w"'^(m^ x) (4.7) 



m2 - D^et 

and 



-^—w^{m\x), (4.8) 







so that 

V„r^(- D^et, a; ) = F(- □.e*, a; ) , (4.9) 

then, by (4.4) and the result in [6], the behaviour of the function (4.8) should already 
be 

x)^- «;^(0, x) log i + 0(1) , e ^ (4.10) 
w^(0, a;) = finite. (4.11) 
The vector current (4.8) is obtained from (3.15) as an expansion: 

x) = x) + x) + 0[K3] (4.12) 

where and Jg are the contributions of the actions S{2) and >S'(3). For the contri- 
bution of the action S{2) one finds by using eq. (3.16): 

I 1 

^reux) = (7l(-nret) + 272(-nret))V'^i? + 2[Va,7l(-°ret)]i?"''+ (4.13) 

+ 2[V^72(-n.et)]i?, 

and 

m.^) = (71(0 + 272(0) V''i?(x) + 0[^^] (4.14) 

whence the constraint (2.14) for the second-order form factors immediately follows. 
Should the constant a in (2. 12), (2. 13) be nonvanishing, the commutator terms in 

22 



(4.13) would contribute to the constraint condition for the third-order form factors. 
Since, however, a = as discussed above, these commutator terms are 0(l/r^) as 
seen from their spectral forms. 

The constraint condition for the third-order form factors is thus of the form 



a,,t,x) + 0[^]=0(r-\x)) , x^I+ (4.15) 

and 

l!;{^,x) = VJ^^i^,x) + 0[^'] (4.16) 

since, in this case, the commutator of with ,^ = — □ contributes to 0[^^] already. 
On the other hand, with the field-theoretic form factors [15] one obtains 

/r(C,^) = ^A^^{x) + {\ogOB''^{x) + 0{e): (4.17) 

I^itx) = -^A^(x) + O(\og0,^^0, (4.18) 

A^{x) = V«A"^(x) +0[3?3] (4.19) 

where the coefficients A^{x), etc. can be expanded in some vector basis second-order 
in the curvature: 

A'^ix) - tr^^p(ni, □2)3?i5R2'^(p) + 0[^'] . (4.20) 
p 

Examples of the basis structures in (4.20) are 

3f?i3?2'^(l) = WRi-R2i, (4.21) 

M2^(2) = VRf ■Va'^pR2i , (4.22) 

^,^^^^{3) = 7^fV^7^2a/3, (4.23) 

3fJi3fJ2'^(4) = Rr^aP2: (4.24) 

etc. where all curvatures of the set (2.2) participate, and the trace in (4.20) refers 
to the matrices i, P, etc. (cf. a similar construction of the basis of invariants in 
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[15]). We do not present the basis in (4.20) in full although it is important to have 
it in full for obtaining the results below. In order that (4.15) hold for any choice of 
the in-state ^ there should be 

A^'ix) + 0[^^] = (4.25) 

and hence 

Ap{n,,a^)=0. (4.26) 

Eqs. (4.26) are the constraints to be satisfied by the coefficients of the l/D 
asymptotic behaviours of the third-order form factors. Let us introduce a notation 
for these coefficients: 

!,(-□!, -Ds, -Da) = ^Flin^, a,) + 0(log(-ni)) , Qi ^ -0 (4.27) 

!,(-□!, -Da, -Dg) = ^F^{D^, Di) + 0(log(-n2)) , Da ^ -0 (4.28) 

'-'2 

r,(-ni, -Da, -Dg) = °2) + 0(log(-n3)) , Os ^ -0 (4.29) 

'-'3 

where there appear functions of two variables FJ^ with m = 1 to 3 and i = 1 to 29. 
The functions Ap in (4.20) are certain linear combinations of the FJ^. By an explicit 
calculation with the action 5(3) in (2.8) and [15] one can work up these combinations 
to see if they vanish. The commutator and potential curvatures appear in the basis 
in (4.20) but the contributions of 57^^,^ and 6P to (3.14) may, in this calculation, 
be omitted since all with acting on TZfj^u or P vanish [12,15]. This makes 
it possible to carry out the check of asymptotic flatness for a generic quantum field 
model. 

^ With the specifications made in sec. 2, this choice boils down to the choice of T^ource 
(2.18). Variations in T^ource induce variations in the curvature of the solution. Eq. (4.15) should, 
therefore, hold for any configuration of the curvature. It is also important for inferring (4.25) that, 
by the construction of the curvature basis in the action [5], /^^(^,x) can have no total derivative 
terms of the form □X"^(^,a;) or \/'^"X^'^^,x). Hence A°'''{x) can have no such terms. 
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The results are as follows. Of 3 x 29 functions only 21 in the table of ref. [15] 
do not vanish and are not related to each other by the symmetries of the form factors. 
With the 21 nonvanishing the expansion (4.20) gives rise to 14 constraints (4.26) 



which, by linearly combining them, can be brought to the following form: 

fU^i,^2) = ^(□i-n2)F4(ni,n2) , (4.30) 

i^io(ni,n2) = __L(n2-n,)2F|3(ni,n2) , (4.31) 

i^28(°i,°2) = 2F|7(ni,n2) + ^(n2-ni)F2yni,n2) , (4.32) 

i^l4(°i,°2) = ^F^^Di, 02) -^(Di- □2)^2^(02, Di), (4.33) 

FU^i, 02) = -\fU^^, 02) + ^(Di + □2)i^28(°i, °2) , (4.34) 

F,\(ni,n2) = --^{u^-u,)F^^{u,,u^) , (4.35) 

i^23(°i,°2) = -2F|2(ni,n2) + (ni-n2)F27(°i,°2)+ (4.36) 

+ + □2)F28(°2, Dl) + ^(^2 - □ODlF^yni, 02) , 

i^9(°i,°2) = ^(□i + n2)(ni-n2)'F2yni,n2)- (4.37) 

- ^(°? + °2)i^l7(°l,°2), 

i^22(°i,°2) = -jF2\(ni,n2) + ^(ni-n2)'F2yni,n2), (4.38) 

Fl{n^,u^) = i(ni-n2)2F2yni,n2) , (4.39) 

i^i6(°i,°2) = (Di -□2)^2^01, 02) , (4.40) 

i^8(°i,°2) = ^(01 + □2)^2^(01, 02), (4.41) 

^i8(°i,°2) = ^F2\(ni,n2), (4.42) 

i^iUni,n2) = -^F2\(ni,n2). (4.43) 



It is now a matter of a direct inspection to check if the F™ calculated in [15] satisfy 
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these constraints. They do! 

Relations (4.30) - (4.43) leave only 7 independent nonvanishing for which 
one can take the functions F21, F25, F2Q, F|g, F2g. With the exception of 

and Fjy, these functions are symmetric in their □ arguments: ^2^5, F|y and F^g are 
symmetric owing to the respective symmetries of the form factors r25, and r28, 
and the symmetry of F21, is a property of the explicit expressions in [15]. 

By expressing all F/" through the 7 independent ones, one can bring the coeffi- 
cient A°''^{x) in (4.17) to the form 

A^"'{x) = -VaV^X['^"ll'^^l(x) + 0[^^] (4.44) 

in which the fulfilment of condition (4.25) is manifest, and the function K'^°"^^{x) 
which appears antisymmetrized in (4.44) is of the following form: 



= 3trF2yni, □2){4V'^V"i?7'' • W^VaRf + 2(^2 - Oi) [v^^Ri ■ V"^i?f ] - 

- 02(02 - a,)g'^'^[Rf . i?2]}i + trF|8(ni, □2){4V'^i?r • V'^i?2?}i + 

+ trF2\(ni, □2){8V^V'^i?f ■ i?2}i + 8 trF2yni, □2){v^V'^i?f • A} + 
+ trF2^5(ni, □2){4i?f • i?2"" + 2g>'''Rf ■ R^}! + 

+ ^trF2\(ni, D2){nr ■ nf) + trF|7(°i, ^2){9>V^'R\^' ■ V,Rf + 

+ 4 g^^V^R^ ■ V^Rf + 4V('^i?i • V'^Rf + (4.45) 

+ 2(ni - □2)^''"[i?i • Rf] + g^'^VV^Ri ■ i?2}i ■ 

The one-loop expressions for the functions F entering (4.45) are given in [15] but 
one may conjecture that the 7 independent structures in (4.45) and (4.44) (5 in the 
case of pure gravity) is the general result independent of models and approximations 
(although this remains to be checked by repeating the analysis above for the general 
form of the action). As will be seen below, expression (4.45) gets considerably 
simplified when inserted in the formula for the news functions. 
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5 The news functions. Contribution of the 



third-order form factors 



A significance of tlie !/□ asymptotic terms in tlie vacuum form factors is that 
they contribute to the energy of the outgoing gravitational waves. This energy is the 
term with the news functions dCi/du, dC^jdu in the Bondi-Sachs equation (3.1). 

Obtaining the news functions requires solving the dynamical equations already. 
However, there is a short cut : eq. (2.4). We may use the fact that the news 
functions appear as a coefficient of the 1/r behaviour of the Riemann tensor at null 
infinity. Indeed, we have (see appendix A) 



1+ 



r du 



where 



C = Ci + iC2 , 

and is a complex null vector tangent to the 2-sphere S: 

(m, Vu) — (m, V^;) = (m, m) = 0, (m, m*) — —2 



(5.1) 



(5.2) 



(5.3) 



with Vw, Vv and S in (3.1)-(3.3), and m* complex conjugate to m. The contribution 
of the outgoing gravitational waves to the mass loss, eq. (3.1), is then 

2 



At: I 



d_ 
du 



and 



C 



(5.4) 



(5.5) 



1+ 
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On the other hand, the Riemann tensor can be calculated with the aid of eq. 
(2.4). Its 1/r behaviour is then obtained as the leading asymptotic behaviour of the 
retarded Green function. Using (2.4b) we have 

92 



1+ 



(5.6) 



where the quadratic terms are to be copied from (2.4b), and the quadratic terms with 
no overall !/□ factor - all terms in the curly brackets in (2.4b) - do not contribute. 

Expression (5.6) can next be simplified as follows. The derivatives in (5.6) appear 
projected either as V^vVi, or as m'^'Vu- In both cases the projected derivatives can 
be commuted with the remaining factors of Vv and m since the components of V 
and Vm in the null tetrad basis are 0(l/r) at X+. The projected derivatives become 
then acting on some scalar X which behaves like 1/r at T+, and in this case 



VvV^X 



O 



x+ 



d 



1+ du 
(see appendix A). In this way one obtains 



= -2 — X + O- 



(5.7) 

(5.8) 



du' 



(5.9) 



x+ 



Eq. (5.9) can now be integrated over u from —00 to a given point of X"*" to obtain 
the news functions. Since, at u — —00, dC/du — 0, there remains to be shown that 
the expression in the curly brackets in (5.9) also vanishes ai u — —00. This is shown 
in appendix B. As a result one obtains 



d_ 
du 



C 



du. i 



2 "^'"^^n 



□ 



□ 



(5.10) 



2{V,^R''){y,^RP^) + 0[R%} 



x+ 
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where for the Ricci tensor one can use the dynamical equations 

There is, of course, a classical gravitational radiation or, more generally, a ra- 
diation induced by the classical source T^^ource ^l- (2-18). The respective contri- 
butionQ to the news functions (call it dC source! du) is obtained by substituting for 
its Einstein value 

Rcl = 8^ (^source ~ 1^ d'^^Tgource^ (S-H) 

in all terms of (5.10) including 0[-R^]. The remaining contributions in (5.10) stand 
for the gravitational radiation induced by the vacuum stress, 

d d d 

C -;r Cgource ~l~ C^ac ■ (5.12) 

OU OU OU 



x+ 



(5.13) 



By (5.10) and the dynamical equations (2.18)0, 

|^C... = -|^{4vrrm,m4(r:,U 0^)} 
and 

Trac = TZi^) + TZi^) + Om (5.14) 
where TZd'^) and TZd^^) are the contributions of the actions S{2) and S{3) in (2.5). 
In (5.10) and (5.13) there appears an expression of the form 

with some tensor X'^'^, and we have used the fact that a contribution of this form 
with 

X^"" = g^'^X (5.16) 

This method can also be appUed to the classical problems in the gravitational radiation. 
Even this contribution is not purely classical since the metric to be used is the solution of the 

expectation- value equations. The same concerns the flux of T^^^^^^ at X+ in eq. (3.1). It would be 

interesting to consider a case where T^^^^^^ does not induce the gravitational waves but T^^^ does. 
Since, at I~ , the flux components of T^^"^ vanish, eq. (5.13) is valid even without the limitation 

implied in (2.4b) (see sec. 2). 



(5.15) 

1+ 
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vanishes by the orthogonahty relation in (5.3). Another property of expression (5.13) 
which will be used below is that a contribution of the form (5.15) with 



X'"' = VX" + VX'' (5.17) 
is of a higher order in the curvature : 

OlX^" X R] . (5.18) 

Indeed, in this case, 

mpm^ ^X^'^ = 2mpm^V''^X'^ + 2m;3m^^([V^ njix'^) . (5.19) 

By (5.7), the first term of this expression is 0(l/r^), and the remaining term contains 
a commutator. 

To lowest order in the curvature, only the contribution of the action S{2) is to 
be considered. By (3.4) and (3.14), 

TZ(2) = ^^{'^''^air + ^''^Jr-9'''^P^alf- (5.20) 

where 



/r = /r(-n,,„x), (5.21) 

and the explicit form of /2(^, x) is given in (3.16). The first two terms in (5.20) are of 
the form (5.17), and, therefore, their contribution to (5.13) is The third term 

in (5.20) is of the form (5.16), and, therefore, its contribution to (5.13) vanishes. In 
the remaining term of (5. 20), the □ operator kills the Green function: 



+0[R^] , (5.22) 
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and this is the reason why at all the asymptotic behaviours of the form factors at 
small □ are relevant to the gravitational waves. Since the second-order form factors 
7i,72 behave hke log(— □) at small □, we have 

I 1 

X ) - O (-) , x^I+ (5.23) 

and the contribution (5.22) vanishes. Thus the vacuum contribution to the news 
functions begins with second order in the curvature: 

^ Cac = 0[R^] . (5.24) 

At second order in the curvature both S'(2) and S{3) contribute to C^ac- We 
arc presently interested in the contribution of >S'(3) which we shall denote C„ac(3). 
Since the action S'(3) is cubic in the curvature, the accuracy in (5.13) is sufficient 
for calculating this contribution : 

|^a„e(3)= -§;;{^^rm,m^^{T:,^^i3) + 0[R'])}\^^ . (5.25) 

To second order in the curvature, T^'^(3) is of the form similar to (5.20) : 

^rc(3) = ^^{'^'^jr+^'^air-g^'^p'^jf- (5.26) 

with 

I 1 

It^It{-Oret,x). (5.27) 

Again, the first two terms in (5.26) are of the form (5.17), and their contribution 
to (5.25) is Again the third term is proportional to the metric, and its con- 

tribution vanishes. Again the □ operator in the remaining term kills the Green 
function: 

I I 

0[R^] (5.28) 



^ C.„e(3) = 1^ { ^ mfsmjfi- □.et,^ )} 
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1+ 



but this time the contribution (5.28) does not vanish because the third-order form 
factors behave hke !/□ at small □. By (4.17), 



n,x)^ -{^y^'^i^) +log(- n)5'^^(^) , (5.29) 

and only the term with !/□ survives in (5.28). With the expression for A'^^(x) given 
in (4.44) the result is 



d d ( r 

C^ac(3) = — I — m„ 



Expression (4.45) for K'"^''^ can now be simplified by using that (i) all terms 
containing the metric with the indices of K^^"^ can be discarded since any contrac- 
tion among mf^m,,Va^f3 in (5.30) results in a vanishing contribution, and (ii) any 
derivative V with the indices of K'^"'^^ can be treated as in integration by parts 
because the respective total derivative contracts with either m or V, and its con- 
tribution vanishes. As a result, K^°"^l^ in (5.30) can be replaced by the following 
expression : 



(5.31) 



+ 4trF|8(ni, 02) [V'Rr ■ VRi^] i + 8trF|7(ni, ^2) [V^'^i^i^^ • V^Rf] 1 + 
+ 4trF2\(ni, 02) [Rf ■ i?r] i + tr{6(n2 - □i)F29(°i, ^2) + 8^2^(01, ^2) - 
- 4F|7(ni, 02)} [V^V^i^f ■ i?2] i + 8trF2yni, Ds) [V^V^i^f • P2] + 



up 
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6 The vacuum news functions in the lowest 



nonvanishing approximation 



To complete the calculation of the vacuum news functions in the lowest nonva- 
nishing approximation we must consider the contribution of the action <S'(2) : 



r\ ^ vac V / r\ 

ou ou 

- (v^li?^')(v^lT_(2)) - 4(v4^fr)(v4^re(2)) + 



(6.1) 



Here R^i is the notation in (5.11), and in the nonlinear terms of cq. (5.10) we omitted 
the contributions T^ac x T^ac = 0{h^) of second order in the Planck constant. As 
shown above, expression (6.1) is of second order in the curvature; the contribution 
of first order in this expression vanishes. 

The expression (5.20) for T^^{2) completed with terms 0[Rl] is of the form 

1 



(6.2) 



where 

2 6^S{2) 



(6.3) 



with 

'^^'^^^^ ^ 2(ip / dxg"\R^,5^,{-U)Rl^ R5^,{-U)r) (6.4) 

is the contribution of the variations of the lowest-order form factors (see eqs. (3.19), 
(3.24)). When inserting expression (6.2) in the linear term of (6.1), all terms in (6.2) 
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proportional to g^^'^ or to the □ operator can be omitted, and for the contribution 
of the first two terms in (6.2) one can use eq. (5.19). 

In the nonhnear terms of (6.1), the expression (6.2) is needed only up to 0[B?]. 
By using (3.16), we obtain 
1 



nrc(2) 



V'^V"(7i(-n) + 272(-n))i? - □7i(-n)i?'^^ - 
- ^^'^'^°(7i(-n) +472(-n))i?] +0[it:^] . 



(6.5) 



This makes it possible to calculate also the Riemann tensor with accuracy 0[B?] : 

R'^^"' = Rf^" - -7i(-n)V[''V<"i?'^]^> + -^f7i(-n) + 272(-n)) x 

X (^^^V^V" - g'^'^V^V^ - ^'^^V'^V" + fif'^"V'^V'^)i? + 0[B?] . (6.6) 

Here we used eq. (2.4a) which defines also 



Rfi'- = ^(2V'^V<"Kf> - 2V"V<"i?f > + 0[R^]) 



(6.7) 



Finally, by combining the results above, the following expression is obtained for 
the contribution dC^ad'^) / du to the news functions: 



d_ 

du 



d ( V 1 



du [in 

+ 4V"V^^i?f^ • Rl^ - 2V"V^3i?rd • <p + Kci ■ R2CI 



(6.8) 



+ 



V^^R^d ■ V''R2cl - '^V'^RZ ■ ^aR2cl + 



+ (7i(-°2) + 272(-n 

+ —R'^li ■ °^2d - IrZi ■ R2ci\ + {AnfU^^) + 0{h') + 0[R']] 

This contribution involves the lowest-order vacuum form factors in the whole range 
of their dependence on the □ argument. 
The total result is 



1^ C = 1^ Csource + ^ C,„c(2) + ^ C,„e(3) + Oih") + 0[R^] 



(6.9) 
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with (9C„ac(2)/(9-u in (6.8) and 9C„ac(3)/(9-u in (5.30). Note that the news functions 
appear squared in the mass-loss formula (3.1), and the vacuum contribution to 
dC/du begins with second order in the curvature, eq. (5.24). Therefore, in the 
absence of a classical radiation, the energy of the vacuum gravitational waves is of 
order 0[h^R'^] which makes this effect difficult to be noticed in perturbation theory. 
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Appendix A. The asymptotically flat metric at null 
infinity 

The general asymptotically flat metric in a chart covering is built by con- 
sidering a congruence of null hypersurfaces u = const, generated by the light rays 
reaching X"*". The generators are labeled by two parameters 9,Lp taking values on 
a 2 - sphere S, and the himinosity distance r is used as a parameter along the 
generators. The metric is then of the form [21] 

ds'^ = -Vdu^ + 2^dudr + r^fab{dx'' - V'du) {dx^ - U^du) (A. 1) 

where 

a, 6 = 1, 2; = 6^ x'^ = (p ^ 
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and 



fabdx^dx^ = \{e^^ ^e^^)de'^ ^{e''-^ -e^-^)smededip^ (A.2) 

+ l(e-2^ + e-2^)sin2^V 



l = (Vii,Vr)<0, (A.3) 



\2 



-1 . (A.4) 



In this metric, (Vu) = is to ensure that the hypersurfaces u — const, are null, 
(Vx", Vu) = is to ensure that the lines u — const., = const, are null geodesies, 
det/ab = sin^ Q is to ensure that r is the luminosity parameter along these geodesies, 
condition (A.3) is to ensure that this parameter is monotonic, and condition (A.4) 
is to choose the retarded time u coincident with the proper time of an observer at 
large and constant r. 

At the limit of X+(r ^ oo, m = const., Q =const., (p =const. ) the metric behaves 
as follows [21]: 

V-I-^^O(^). (A.=) 
7 + * Ci 



^ = ^-(;^)' 

= — + O (pj) , (A.8) 

* = +0(;^) (A.!>) 

where tM, Ci, C2, A^", B arc functions of (/?, u. 

The Ci,C2 differentiated with respect to u are the Bondi-Sachs news functions 
[20, 21]. In the gauge (A.l) they stand for the radiation degrees of freedom of the 
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gravitational field. The Bondi mass M[u) is obtained by averaging the coefficient 
A4 in (A. 5) over the unit 2-sphere : 

2-7r 7r 



M{u) 



-I 



d'SM 



dip / d9sm9M{9,ip,u) 



47r J An 



Its limiting value at m — > — oo is the ADM mass, and the difference 

dM. 



(A. 10) 



M(-oo) - M{u) = J du(- 



du^ 



(A.ll) 



is the energy radiated away through X+ by the instant u of retarded time ( see [22] 
and references therein). 

Some components of the Ricmann and Ricci tensors calculated in the metric 
(A.l) - (A. 9) at X+ are as follows 

^2 



Rr0 

R. 
R, 



— r 



1+ 



1+ 



x+ 



1+ 



x+ 



x+ 



Ci + 0{1) , 



rsin^-^C2 + 0(l) , 
ou 



rsin^^^Ci + Ofl) , 
ou 



ll + AB 



O 



2iV^ + 



1 



d 



(Ci sin^ e 



d / a 



dip ^sin^ 



(A.12) 
(A.13) 
(A.14) 
(A.15) 
+ 0(1), (A.16) 



2N^sin^e+ J- A(cr2sin2^)-^Cil+0(l) , (A.17) 



1 A 
du 



(Cf + C| + 4i3)+0(l), 



(A.18) 



"'^^By our calculations, the equation Rrr = {2/r){'^~^dr^ — (drj)"^) ( in the present notation) 
given in [21] for the axisymmetric case is in error. We obtain Rrr = (2/r)^'^^9rVE' — 2(9r7)^. It 
follows from the former equation that, for i?^^ of a compact spatial support, ^' = — l + 0(l/r'^). 
With our result, \E' = — l + 0(l/r^). The significance of the behaviour of \1/ is seen from eqs. (A.15), 
(A.18). 
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Rue 
R 
R 



Uif 



1+ 



1+ 



1+ 



2N^ 



2N^ sin^ ^ + J:^ 4 (^2 sin^ 9) - C^\ + O (1 



1 ^ 

1 _a_ 

r du L"^' ""^ " ' sin^ 89^^^""^ dcp 

-du ~'~ ^du ~^ ^du ^) ~'~ sin^ ' (^i' 



sin^ 09 



,(A.20) 

+ 
(A.21) 



When referring to tensors at X+ we always mean their projections on the null 
tetrad Vu, Vf , m, m* introduced in (3.2), (3.3) and (5.3) where v and rria are asymp- 
totically of the form 

1^ 



2r + u + 



r{VJ + ism9Va^) + 0{l) 



(A.22) 
(A.23) 



The null-tetrad components of physical quantities are regular at T"*" i.e. are either 
finite or decreasing like inverse powers of r. This is true specifically of tensors ob- 
tained by the action of the retarded form factors (see below). The null-tetrad vectors 
may be regarded as covariantly constant at X+ since the null-tetrad components of 
their derivatives are 0(l/r). Thus, up to curvature terms. 



— (m^m; + m*ma) + 0[R„ 



(A.24) 



x+ 2r 

and the curvature terms are 0{l/r). Eqs. (5.7), (5.8) are obtained by calculating 
the derivatives projected on the null tetrad in terms of the Bondi-Sachs coordinates. 
Specifically, eq. (5.7) owes to the fact that = 0{l/r) as seen from (A.23) and 

(A.1). 

It follows from the asymptotic expressions above that the null-tetrad components 
of the curvature tensor decrease at X+ like 1/r or faster. Specifically, 

x+ r ^du^ ' du 



(A.25) 
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which is eq. (5.1). The energy flux component of the Ricci tensor at X"*" 



1+ 



(A.26) 



is given in eq. (A. 21). By averaging (A. 21) over the unit 2-sphere one obtains the 
relation 



dM{u) 
du 



x+ 



which, after using the dynamical equations 



1 

2^ 



(A.27) 



(A.28) 



becomes the conservation law (3.1). Here T^^^^l is the total energy-momentum tensor 
which in eq. (2.18) is 

(A.29) 



rp^lV rpjXV I rpjXV 

total source ' vac " 



Appendix B. The retarded Green function in the 
past of X+ 

To lowest order in the curvature, the retarded operator !/□ acting on an arbi- 
trary tensor source X"^ ' "* is of the following form: 

= ^ / dx-g'lH{a{x,x))gX{x,x) . . . gX{x,x)X^'-^^{x) + 0[X x 3fJ] 

post of X 

where (j{x^x) is the world function [23], or geodetic interval biscalar [24], g'^^{x,x) 
is the geodetic parallel displacement bivector [24], and the integration point x is in 
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the past of the observation point x. Here and below, the bar over a symbol means 
that this symbol refers to the point x. 

It follows from a comparison of the equations defining 51"^ (a;, x) with the ones 
defining (t{x,x) that, up to the curvature terms, the parallel displacement bivector 
can be calculated as follows: 

g\{x, x) = -V"V,,a(a:, x) + 0[R] (B.2) 

whence it also follows that 

V,g\{x,x) = 0[R]. (B.3) 

If with i = 1 to 4 are the vectors of the null tetrad, then the null-tetrad compo- 
nents of the tensor (B.l) are obtained by calculating the contractions 

lt{x)g^-,{x,x) . (B.4) 

By using eq. (B.2) and a perturbative expression for the world function, it is easy to 
see that, when x tends to X+, and a; is in a compact domain, the contractions (B.4) 
remain finite. Hence the null-tetrad components of the tensor (B.l) decrease at X+ 
like 0(l/r) - the fact assumed in the main text. The expression for the massive 
retarded Green function is similar to (B.l) [6]. Therefore, generally, the null-tetrad 
components of tensors obtained by the action of the retarded form factors are regular 
at I+. 

As seen from (B.l), we are always dealing with some scalar source 

y{x) = gXi^,x) . ..gX{x,x)X'^-^'^{x) (B.5) 

which may depend parametrically on the observation point but it suffices to consider 
the action of the Green function on a scalar : 

- ^X"i-"'=(x) = ^ j dxf^^6{a{x, x))y{x) + 0[X x 3fJ]. (B.6) 

past of X 
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The integration over the hght cone in (B.6) includes subintegrations along the light 
rays coming from X~ to the observation point x : 



(B.7) 



Here L is a generator of the past light cone of p is the luminosity parameter along 
L, and the measure in (B.7) is asymptotically of the form 



d\x{p) 



p—KX> 



dp ■ p 



We shall, therefore, assume that the source decreases at I like 



X 



ai...«fe 



o 



(B.8) 



(B.9) 



Another important assumption [6] is analyticity of the source in time including 
the past timelike infinity The real sources appearing in the expectation- value 
equations are built out of the curvature, and the condition of analyticity implies 
in particular that, at i~, the metric becomes asymptotically static. This should be 
provided by imposing the respective requirement on T^^^^. By analyticity, the limit 
r ^ oo of the source at coincides with its limit in the past of X" : 



where v is the advanced time along I . Hence, by (B.9), 



(B.IO) 



:b.ii) 



For the calculation of the integral (B.6) at x — > X+ we may use the Bondi-Sachs 
frame (A.l). For the past of this is safe even if the metric has closed apparent 
horizons [2,3] since no one of these will be encountered by the hght rays emitted 
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sufficiently early. To lowest order in the curvature, the world function is then of the 
form 



u — u. 



a{x, x) = —{u — u){r — r -\ — ) + rr(l — cosa;) + 



(B.12) 



where is the world function on the 2-sphere : 



cos u; — cos 9 cos ^ + sin ^ sin 6 cos(</7 — </?) . 



(B.13) 



By solving the equation a-{x, x) = with respect to u and choosing the solution 
which corresponds to the past light cone of x,u — f, we obtain 



□ 



-X 



1 /" — /" — — 0(7 —1 — 

, e,^,u)^ — j (fs J dff^ (Vm, Vf) • — y{f, e, ^, u) 



+ 



u=f 

+ 0[Xx^ (B.14) 



and 



where 



1+ 



^Q(9,ip,u) + 0(\ 



(B.15) 



Q{e,ip,u)^-— I d's I dfiiogf)—[f'yif,e,ip,u)l^^^]+ o[Xx^ , (B.ie) 



f* ^ u - f {1 - cosuj) + 0[R] , 



(B.17) 



and the equation u — f* is the equation of the limiting light cone of the point x 
at T+. Here we wrote = f^dlogf/df and integrated by parts for being able to 
consider sources decreasing at I~ like At the limit w — > — oo, the source in 

(B.16) turns out to be at — > — oo. By the assumption of analyticity, we then have 



Q{0,(p,u) 



Qo + -Qi + 

u 



(B.18) 
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with 



1 

An 



d r - 1 

d'S / rff(logf)— f33^(f, e, -oo) + 0[X X 3fjl 



(B.19) 



etc., and the convergence of the integral in (B.19) is now owing to (B.ll). 

Eqs. (B.15)-(B.19) make manifest the fact which is more general than the ap- 
proximations made. Namely, the integral (B.l) in the past of T"*" involves only the 
source at i~ . Under the assumption of analyticity at i~, this integral be- 

comes a static Coulomb potential. One can show that it preserves this form also at 
two other limits : at the past null infinity and spatial infinity (i°). 



□ 



X~ or orX+,u- 



^(l + O), 
r 



;B.20) 



Consider now equation (5.9) for the news functions which is valid under the 
assumption made in sec. 2 that the flux components of the Ricci tensor at X~ 
vanish. Since, in this case, R^^ at X~ is 0(l/r^), the integral implied in (l/D)/?^'" 
converges. Moreover, the derivatives of this integral in the null-tetrad basis behave 
at X~ like 



(v.-- 



O 



:b.2i) 



(see below). Therefore, the nonlinear additions to R^^ in eq. (5.10) are 0(l/r^) at 
. We then have 



mprriy — 



□ 



□ 



□ 



□ 



ig(^,^,..) + 0(l) . (B.22) 



Since, at Z+, the flux components of i?^'' do not vanish, the integral {l/U)R^^'^ behaves at T+ 
Uke (log r)/r. The (log r)/r asymptotic terms vanish in expression (2.4b) owing to the presence of 
the antisymmetrized derivatives, and in expressions (5.10), (5.13), (B.22) owing to the contraction 
with mfjm^. 
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We need to calculate the limit 



d 

lim — q{9,^,u) 

u-»-oo (ju 



(B.23) 



which serves as an initial datum for eq. (5.9). By analyticity of the metric at i , 



= ?o + - gi + 



(B.24) 



and the limit (B.23) vanishes owing to the presence of the derivative d/du. 

Eq. (B.20) can also be used to prove that the solution of the Bianchi identities 
(2.3) with zero initial data for the gravitational field at X~ is expressed indeed in 
terms of the retarded Green function. This is easily seen in expression (2.4b). Since, 
by (B.20), 



-i(i?^^ + 0[i?^]) 



+ 



(B.25) 



and Qq^ does not depend on the advanced time along I~ , the derivatives in (2.4b) 
enhance the power of 1/r. Eq. (B.21) is valid for the same reason. 

The proof is more involved in the general case, where the fluxes of R'^'^ through 
I~ are nonvanishing, since one has to address eq. (2.4a). In terms of the null tetrad 
at 2" 



y[/iy<a^i/]/3> 



^-^(V,K)(V.u)VW'^i?^l/^> + 0[^). (B.26) 



Since the flux components of R"^ at X~ are proportional either to Vv or to V^-y, 



they cancel in (B.26) owing to the antisymmetrizations : 



y[/^y<a^i/]/3; 



o 



I- 

Thus the source in (2.4a) satisfies condition (B.9), and we have 



(B.27) 



(B.28) 
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with Q%^^^ constant along X . By (B.19) and (B.5), 



1 



and, up to is a total derivative: 



y 



;b.29) 

(B.30) 

(B.31) 
fB.32) 



where use is made of eq. (B.3). Since, in this total derivative, the time derivative 
vanishes owing to the asymptotic stationarity of the metric at and the angle 
derivatives vanish in the integral / d^<S, we obtain 



-iz I (i'S I df{\ogf)^ 



dr 



d 



+ 0[Rl 



(B.33) 



Since the power of decrease oi Z'^ at r ^ oo is at least (9(l/r^), we may integrate 

i 

by parts in (B.33) to remove the integration over f completely : 

Qf^'' = ^lim i- j ci25(4^[''V''^'"V^^V^f) (f2V^^,^[_) + 0[R^,] . (B.34) 

Only the terms in Rj^p that decrease hke can survive in (B.34). However, the 
integrand in (B.34) contains one more derivative of ^p^, and, again, the respective 
time derivative vanishes at i~ . Therefore, this integrand is 0{l/f) (actually 0(l/f^), 
see below), and the hmit (B.34) vanishes. The presence of an extra derivative is in 
fact not essential. By (B.IO), the sequence of limits in (B.34) can be replaced by the 
limit of the same quantity in the past of X~ . It is then seen that the constant (B.34) 
hangs solely on the flux components of the Ricci tensor in the past ofX^. However, 
the metric in the past becomes asymptotically static, and for a static metric the 
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fluxes of Ryp tlirougli X are absent as seen from tlie counterparts at X of eqs. 
(A. 18) - (A.21). We have, therefore, 



{^i^P ,_) 



and 



Hence 



Qf"'' = . 



R 



I- 



;b.35) 



(B.36) 



;B.37) 



which proves that, in the retarded solution, there are no incoming gravitational 
waves. 
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